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INTRODUCTION 


The study of airframe and engine-inlet integration has recently become even more 
important in the light of current fighter and supercruiser aircraft design. These 
configurations are characterized by fuselages and canopies with distinctly nonaxisym- 
metric cross-sectional geometries and thin, low-aspect-ratio, highly tapered swept 
wings. Inlet external surfaces are often well integrated with the fuselage (ref. 1), 
and ^the inlets may be located almost anywhere on the fuselage (ref. 2). Such air- 
craft configurations may also include, in more complicated cases, forward -mounted 
canards, wing strakes, leading-edge extensions (LEX's), fins, and stores hung in a 
multitude of combinations. 

Inlet performance is determined to a great extent by the flow field just ahead 
of the inlet entrance plane. Satisfactory inlet performance usually requires that 
flows at the inlet face be high in total-pressure recovery and that this recovery and 
the local flow inclination be relatively insensitive to aircraft flight conditions. 

The effects of these configuration parameters on the flow into the inlet are 
many. The inlet is subject to a number of disturbances from the aircraft forebody 
alone, including large flow variations, substantial boundary- layer growth, vortex 
shedding, and shock impingement. If the inlet is shielded by the wing of the air- 
craft, further complications may arise through the complex interaction between the 
wing and forebody, the effects of which are highly three-dimensional. In such cases 
it is usually possible to adjust the inlet location and the wing position to optimize 
inlet performance, but not without sacrifices. For example, burying the inlet well 
behind the leading edge of the wing near the fuselage will almost certainly improve 
the inlet performance at high angles of attack, but usually at the expense of per- 
formance at substantial angles of sideslip (ref. 3). 

Other contributors to aircraft inlet flow problems are the various protuberances 
that are placed on the aircraft fuselage or wings. Canopies, LEX*s, and canards, in 
their traditional positions, cause the most serious problems for engine inlets at 
high angles of attack or sideslip or both (ref. 4). Low- or side-mounted inlets are 
more likely to be disturbed by stores hung from the fuselage or wing or by ventral 
fins, although the latter are usually placed far enough aft to be no problem. 

The purpose of this study was to assess the capability of two numerical flow 
prediction methods to predict forebody flow fields near the engine inlet on realistic 
fighters. For this study only inviscid prediction methods were considered. The 
target of this study was to determine the shortcomings of these methods, vrtiich could 
not be expected to predict viscous effects. The fighter speed range from Mach 0.9 
to 2.5 and angles of attack from 0° to 25® were of interest, and of prime importance 
was the realistic modelling of the aircraft geometry. The ease of use of the com- 
puter programs that implemented each method and their economy of operation was also 
considered. 

Predictions of aircraft forebody flows are valuable not only to the designer but 
also to the wind tunnel test engineer (ref, 5). The results are presented in a form 
particularly useful to both, namely, flow-field contours, at the hypothetical inlet 
entrance plane, of local angle of attack, angle of sideslip, and Mach number. 



Two other subjects, which will not be included in this investigation, must be 
mentioned in the interest of completeness: cyclic flow distortion and spillage. The 

first phenomenon may have a strong adverse effect on inlet pressure recovery, partic- 
ularly during maneuvering flight (ref. 6), and is extremely configuration dependent. 
Inlet spillage, because of its deformation of the aircraft flow field, can seriously 
affect the aerodynamic efficiency of the aircraft. Recent theoretical predictions of 
this phenomenon have agreed well with wind tunnel data (ref. 7). 

SYMBOLS 

Cp specific heat at constant pressure 

c^ specific heat at constant volume 

M Mach number 

P In p 

p pressure (normalized with respect to p^) 

S entropy (normalized with respect to c ) 

V , 

T temperature (normalized with respect to T^) 

u,v,w velocities correspo nding to the x, y, and z directions (normalized 
with respect to \jp^/p^) 

x,y,z Cartesian coordinates in the axial, transverse, and normal directions 

(normalized with respect to an arbitrary length) 

a angle of attack, tan"'^(w/u) 

3 angle of sideslip, tan“^ (v/u) 

Y ratio of specific heats, c^/c^ 

p density (normalized with respect to p^) 

()) perturbation potential function 

Subscripts : 

i local conditions 

00 free-stream conditions 

x,y,z partial derivatives with respect to the appropriate direction or directions 

A bar over a symbol indicates a dimensional value. 

NUMERICAL PREDICTION METHODS 

Because of the differences in the simplified equations of motion for supersonic 
and transonic flows, almost all the numerical methods investigated were specifically 
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for one speed regime or the other. During the course of this study, certain criteria 
were used to assess the applicability of each of the computer programs. Solution 
accuracy, integrity of the geometric model, operational ease, and running expense 
were the most important criteria. In the final analysis, for supersonic cases, a 
three-dimensional Euler equation marching code, STEIN, was chosen over a surface 
paneling code, the PAN AIR Pilot Code. Although the latter program was applicable 
to a wide range of geometric configurations and had reached a high level of technical 
development (ref. 8), recent studies (ref. 9) have shown that it is less accurate 
relative to the STEIN code and is substantially more expensive to operate. 

A similar situation developed in the search for a transonic code. Ultimately, a 
small-disturbance code, WIBCO, was chosen over a full-potential code, FLO-30, which 
is the latest member of a family of FLO-codes (refs. 10 to 12), each of which is 
capable of handling a more complex wing-body configuration. Investigation of the 
most recent FLO-30 calculations (ref. 13), however, indicated that the method by 
itself was capable of solutions for fuselages of only moderate complexity, particu- 
larly if the wing was of low aspect ratio. Since many of the fighter cross sections 
varied considerably from an axisymmetric shape, it was decided that WIBCO would be 
the better choice, in spite of the small-disturbance approximations applied in the 
calculations . 

Although neither WIBCO nor STEIN had been developed for application to fighter 
aircraft configurations, it was anticipated that modifications to the codes as might 
be necessary could be efficiently carried out because of their usage in many areas 
and the availability of complete documentation. 

One other factor weighed heavily in the choice of the STEIN and WIBCO codes: 
the similarity of the aircraft model geometry input. Both codes have slightly dif- 
ferent versions of the QUICK-geometry methodology developed by Vachris and Yaeger 
(ref. 14). In this method the aircraft cross sections are described at appropriate 
stations by specifying control points at each cross section and the type of curves to 
be used to connect the points. After the details of each cross-sectional geometric 
model are specified, the aircraft body lines are described in a manner similar to the 
cross sections, piece by piece, with lines, ellipses, parabolas, and so forth. The 
result is an analytical body model that allows a quick calculation of surface points 
and slopes over the entire length of the body. This method requires the body coordi- 
nates to be single-valued polar coordinates. The wing, if present, is considered as 
an outgrowth of the body for the STEIN input, whereas the wing is a separate specifi- 
cation of chord sections at different span locations for the WIBCO input. 

The task of describing the aircraft geometry by this method can become time 
consuming. It was eased substantially for the present study through the availability 
of an interactive graphics system developed by Adams (ref. 15). Examples of the use 
of this system in defining a typical fuselage cross section and a body line are shown 
in figure 1(a). The resulting geometric model, plotted as a series of cross sec- 
tions, is shown in figure 1(b). 


Supersonic Code - STEIN 

According to reference 16, STEIN was originally formulated to provide inviscid 
solutions for flows in supersonic or hypersonic conditions for a wide variety of 
vehicle configurations throughout a large angle-of-attack range. All forms of small- 
perturbation techniques, valid only at low Mach numbers and low angles of attack, and 
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Newtonian theory, valid only at high Mach numbers, were rejected as being inadequate 
for the problems at hand. Thus, STEIN was developed as a finite-difference marching 
solution of the three-dimensional Euler equations with shock-fitting techniques 
applied to selected types of shocks. Marching solutions carry an implicit require- 
ment that the flow throughout the field of interest be supersonic. As will be seen, 
this requirement affects the range of Mach numbers and the type of configurations for 
which the STEIN code is applicable. 

The geometric model that is input to the code is in the QUICK-geometry format. 
Wings, canards, and other protuberances are defined simply as outgrowths from the 
fuselage shape at each cross section. The input that the code requires as well as a 
summary of the QUICK-geometry system are presented in reference 17. 

The solution "marches” along an axis of the vehicle from a starting solution 
near the nose toward the tail of the aircraft. Two options are allowed for the 
starting solution, one for a blunt nose and one for a pointed nose. The fighter 
aircraft in this study were configured with pointed noses, and thus the second option 
was used. A conical supersonic flow is assumed up to a certain fuselage station, the 
"starting plane," beyond which the marching solution takes over. Downstream from the 
starting plane, the region of computation is bounded by the body and by the bow shock 
of the vehicle. This area is first transformed, conformally, to a mapped space in 
which the geometric cross sections become "near circles.” The computational space is 
formed from the mapped space by normalizing the radial and circumferential distances 
between the body and the bow shock, with the portions of the plane of symmetry above 
and below the aircraft becoming the upper and lower limits of the computational mesh. 
The normalization results in a rectangular computational grid which transforms into a 
grid in the physical space that is finely spaced in regions of highly convex body 
curvature • 

Shocks that form within the computational grid are considered to be of two 
types: "cross-flow" shocks, which are predominantly radial, and "wing" or "canopy” 
shocks, which are essentially at a constant radius in the mapped plane. Figure 2 
shows a typical shock configuration generated by a winged aircraft in which there 
are bow, canopy, and wing shocks. All these imbedded shocks are located by monitor- 
ing local pressure throughout the flow and determining the point of maximum gradient 
based on a cubic polynomial curve fit. Locations of inverse pressure gradients 
arbitrarily near zero are designated as shock points. The transformation from the 
mapped to the computational space is then adjusted so that the shocks become mesh 
lines, with the Rankine-Hugoniot relation being applied across each shock. As the 
code progresses downstream, each shock wave is followed until its demise or until it 
merges with another shock. A typical representation of the mesh as it appears in 
the physical plane is shown in figure 3. Note the manner in which the canopy and 
wing shocks are fitted into the mesh between the body and bow shock, ensuring that 
the shock locations are single valued throughout all the planes. 

In regions of continuous flow, the Euler equations in the physical plane in a 
marching form are 


UP + YU = -(vP + w P + Yv + Y^ ) 
XX y z y z 


( 1 ) 


UV = -(w + wv + TP ) 

X y Z y 


( 2 ) 
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( 3 ) 


uw = 

-( vw 

+ ww 

X 

y 

z 

TP + 

uu = 

-(vu 

X 

X 

y 

uS = 

-(vS 

+ wS ) 


X X z 

with a corresponding ideal-gas equation of state 


(4) 

(5) 


In T = P + — S (6) 

Y Y 

Upon transformation to the computational plane, the derivatives of the five dependent 
variables P, u, v, w, and S appearing on the right-hand sides of the equations 
are evaluated explicitly, and the solution is stepped by using a MacCormack two-level 
predictor-corrector scheme (ref. 16). This procedure is varied at body and shock 
points. Use of the MacCormack scheme and strict control of the step size ensures 
second-order accuracy in the marching direction for a linear system of equations, and 
truncation error is assumed correspondingly small for the present system. 

The version of the STEIN code used herein includes the most recent modifications 
to the capabilities of STEIN (ref. 18): inlet mass ingestion, subsonic axial Mach 
number, improved conformal mappings, and sideslip. Only the second and third modif- 
ications were used in this study. 

Transonic Code - NIBCO 

The WIBCO transonic code was developed by Boppe (ref. 19) primarily to apply the 
solution of the small-disturbance potential equation to arbitrary wing and body geom- 
etries. Recognizing the increasing complexity of traditional grid transformations as 
configurations become more three-dimensional, Boppe avoids these problems by imbed- 
ding fine Cartesian grids into an overall coarse grid in regions where more flow 
detail is required. The concept is shown schematically in figure 4 for wing and body 
fine grids imbedded in the global crude grid. The crude grid in the physical space 
is stretched in all directions to infinity (except at the plane of symmetry) accord- 
ing to the method of reference 20. Stretching in the z-direction is also propor- 
tional to the tangent function, whereas in the y-direction the function of choice is 
the hyperbolic tangent. 

The wing and body fine-grid systems are constructed to totally encapsulate their 
portions of the geometry and to provide computations oyer a much smaller area of the 
flow. These two fine-grid systems overlap and transfer information to each other, as 
well as to the crude-grid system, during the course of the iterations. It should be 
noted that although the body fine-grid system is a regular Cartesian grid, the wing 
fine-grid system is swept and tapered according to the planform shape. 

The geometry system used for WIBCO is the same as that used for the STEIN super- 
sonic code, with two differences: first, WIBCO uses QUICK-geometry only for the body 

and with several options removed; and second, the wing is defined separately as wing 
sections at different span stations. Once these two limitations are considered, 
however, a careful design of the STEIN geometry can be transferred with only minor 
changes to be used as WIBCO input. 
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The small-disturbance equation used in WIBCO, including three terms that are 
usually present only in the full-potential representation, is 

fi - M ^ - (Y + 1 )M ^(j) - ^ ^ M ^(f> ^l(fi - 2M ^(}) 4> 

00 00 2 00 J XX 00 ^y^xy 

+ [l - (Y - ](j> + <|) = 0 (7) 

L 00 yy zz 

The 4>y<j>xy ^x^yy terms are included to facilitate the resolution of shock 

waves on wings of large sweep angles, and the 4>x^4^xx term is included to provide a 
better indication of when the equation changes type from elliptic to hyperbolic and 
vice versa. 

The finite-difference approximations are straightforward. Central differencing 
is used throughout except in areas of local supersonic flow, in which upwind differ- 
encing is used for most of the second derivative terms. In keeping with the near- 
i sen tropic nature of the flow, nonconservative difference operators are used, 
although it is acknowledged that results will become less accurate with increasing 
shock strength. For a wing-body configuration, the solution begins with an arbitrary 
number (typically 100) of successive line-overrelaxation sweeps of the crude grid to 
provide a starting solution for the fine-grid systems. The second phase of the solu- 
tion involves a sweep of the wing fine-grid, the body fine-grid, and the crude-grid 
system with appropriate updating of overlapping areas. Approximately 80 second-phase 
iterations are usually required. Since none of the grid systems are body or wing 
fitted, boundary conditions are applied at mesh points nearest the actual surfaces. 
Corrections are applied at these points for wing-surface slope and body displacement 
as well as for local flow inclination. 

The version of WIBCO used herein is the basic wing-body code. Further capabili- 
ties have been added (ref. 21), including the fine-grid system applied to pylons, 
nacelles, and winglets, as well as a scheme for modelling inlet spillage and exhaust 
interference effects. 


EXPERIMENTAL DATA 

One of the problems in obtaining good inlet flow-field data for realistic air- 
craft forebodies is the sheer size of the models necessary to provide sufficient 
geometric detail. The problem is made much worse if flow-field studies are extended 
to high angles of attack and sideslip. At least one imaginative method has been 
tested (ref. 22) to provide realistic inclined flow to an inlet while eliminating the 
need for a massive forebody model in the wind tunnel. On the other hand, a great 
amount of data on large models was obtained in the early 1970*s in Project Tailor- 
Mate, a study of various aircraft forebody shapes and their effects on engine-inlet 
flows. The first part of the project, which was specifically concerned with 
forebody-alone flow fields, provided the data used for comparisons herein. These 
forebody tests were performed without engine-inlet simulators on wing-body combina- 
tions representative of three types of fighter aircraft inlet configurations: side- 

mounted, wing-shielded, and body-shielded. 

The first two of these configurations are considered in this paper, and sketches 
of the full aircraft configurations are shown in figure 5. It was felt the third 
configuration would not be as challenging as the first two. The supersonic and 
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transonic wind tunnel tests were performed in the 16S and 16T wind tunnels at AEDC 
(refs. 23 to 25). The actual tests used in this report were run without engine 
nacelles. Three types of data were taken during the forebody tests: cone-probe 
flow-field data, body pressure data, and boundary-layer data. The predictions in 
this paper are concerned only with the forebody flow field near the inlet. Thus, 
comparisons are made with the first type of data and are presented in a format offer- 
ing direct comparison with the data in references 23 to 25. 

The tests were run at Mach numbers of 0.9 to 2.50 and at angles of attack of -5® 
to 25®. Although the quality of the data is considered to be very good, post-test 
examination of the cone-probe data (ref. 26) has indicated that one of the five cone 
probes was faulty. On examination of the original tabulated data, it was decided to 
omit data from this probe. Surber and Sedlock (ref. 26) estimate the data from the 
other probes to be accurate to ±1®. One other problem was the lack of detailed draw- 
ings of the aircraft configurations. As a result, the numerical models generated for 
this study were based on rather small scale cross sections presented in references 23 
to 25. However, the inaccuracies thus introduced appear to be no greater than those 
introduced by the approximations made in defining the aircraft geometry analytically, 
as required by the STEIN and WIBCO codes. 


RESULTS AND DISCUSSION 

A partial summary of the results herein is presented in reference 27. The 
results of this paper are summarized in table 1 . 


Side-Mounted Inlet Configuration 

The fuselage geometry generated with the STEIN code for the side-mounted inlet 
configuration (fig. 5) is shown in figure 6. The density of the circumferential grid 
was increased about halfway down the fuselage, as was the density of the radial grid 
(not shown in the figure), in order to provide a more detailed mesh in the region of 
the inlet. The WIBCO model, being specified by the same QUICK-geometry as the STEIN 
model, was very similar to the model in figure 6, with the density of points in the 
body-fine grid system tending to be higher than those of the STEIN model over the 
entire body. 

Representative cases for the STEIN code were chosen to be those at Moo = 2.5. 
Solutions were successful at this Mach number for angles of attack up to 15®, as 
shown in figures 7 to 10. These figures show comparisons of predicted and experi- 
mental contours of local angle of attack, local angle of sideslip, and local Mach 
number. Beyond a = 15®, the inclination at the nose of the aircraft was too high 
for the conical flow starting solution to be successful, indicating that the actual 
flow contained a detached bow shock. It is conceivable that a blunt-nose starting 
solution would be capable of overcoming this difficulty, but that approach was not 
attempted in this study. 

Examination of figures 7 to 10 shows that both the predictions and the experi- 
mental data behaved consistently as the angle of attack was increased from 0® to 
15®. The agreement of the prediction and the experimental data of the flow inclina- 
tions is good, whereas the agreement between Mach numbers is consistently poor. The 
theoretical local angle of attack is consistently 1 ® high, and the theoretical local 
angle of sideslip is consistently 1 ® low (inboard) , perhaps reflecting the estimated 
±1 ® probe accuracy stated previously. 
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When the Mach number was lowered, the STEIN code was capable of producing a 
good flow solution at Moo = 2.2 and a =5 0® and (fig. 11), which compared as 
favorably with experimental data as the M«j = 2.5 solution (fig. 7). No solution 
at Moo = 1 *6, even at a = 0®, was possible, however, as the configurations and 
flight conditions once again produced a detached bow shock. 

Similar good results were obtained with the WIBCO code in the transonic case, 
shown in figures 12 to 17. In general, local angles of attack and sideslip agreed 
to within a couple of degrees, with discrepancies growing as the angle of attack 
increased to a = 25®. At this highest angle of attack (see fig. 17(b)), the 3^ 
contours show that in the experiment the flow showed a much stronger outflow on the 
lower (windward) side of the body and a much weaker inflow on the upper (leeward) 
side of the body than the prediction, which can be taken as an indication of viscous 
effects and possible flow separation. This can be seen in the Mach number contours 
also (see fig. 17(c)), in which the experimental flows did not accelerate around the 
upper shoulder of the model as the inviscid prediction of the WIBCO code shows. 

The WIBCO code was able to produce solutions for this configuration at low 
supersonic Mach numbers of 1.2, 1.4, and 1.6, and figure 18 shows the results at 
Moo = 1.6. Local flow angles show good agreement. A WIBCO solution at M» = 2.2 
(fig. 19), well above what is usually considered transonic flow, shows very good 
agreement with the experimental data and the STEIN results (fig. 11). 


Wing-Shielded Inlet Configuration 

Figure 20 shows the geometric model for the STEIN code for the wing-shielded 
inlet configuration. As before, the coordinate mesh density was increased about 
halfway down the model. The geometric description shown is one composed entirely of 
cross-section descriptions up to the area of interest. For the transonic cases, how- 
ever, the overall sensitivity of the flow to a much longer portion of the aircraft 
geometry dictated a different geometric description. The aircraft model used for the 
WIBCO code was much longer, and the wing is described in its entirety as a spanwise 
sequence of cross sections. Attempts made to model the geometry in the same manner 
as for the supersonic case proved to be unsatisfactory. 

The supersonic results of the STEIN code were obtainable only at a = 0® and 5® 
at Moo = 2.5 because of large subsonic areas near the wing-fuselage juncture. The 
results are shown as figures 21 and 22. 

Shielding the inlet under the wing provides a substantial lessening of the 
effects of the local angle of attack, as can be seen from comparing figures 8(a) 
and 22(a). Although quantitative agreement of the inclinations does not seem to be 
very good, it is significant to note that angular variations are small, and minute 
changes in local flow angles can cause substantial changes in contour placement. 

Wind tunnel data obtained at a fuselage station slightly upstream of the station 
shown in figure 21 (a) agree much better with the predicted contours even though the 
magnitudes of the angles involved are essentially the same. At angles of attack of 

10® and above the STEIN code encountered difficulties. The problems occurred at the 

sharp leading edge of the wing, where a local two-dimensional flow solution is used 
for the shock wave. The flight conditions and the geometry of the wing called for a 
detached shock, a situation that was unattainable with the present algorithm. Modi- 
fications to the previously described shock-fitting scheme could possibly remedy this 
situation, as could the redefinition of the leading edge of the wing as blunt. The 

results at a = 5® are reasonable, with discrepancies near the body attributable 
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once again to inaccuracies in the geometric description and lack of viscous capabil- 
ity in the prediction code. No solutions were possible for this configuration for 

Mach numbers of 2.2 or below. 

The WIBCO code, run for cases at M^o = 0.9, was more successful at the higher 
angles of attack. Results presented in figures 23 to 26, for angles of attack from 
0® to 15®, show the flow patterns developing consistently. Qualitative agreement 
with the data is reasonable, with flow inclination errors on the order of 3® maximum 
at the highest angle of attack. As in previous cases shown for both configurations, 
both geometric modelling and viscous effects are responsible for discrepancies in the 
results. There is evidence of a large boundary layer moving from the fuselage to the 
lowe^r surface of the wing, especially at the higher angles of attack. 

Results of the WIBCO code at a = 0® for three higher Mach numbers (up to 

Moo = 1-6) are shown in figures 27 to 29. These figures and figure 23 show agreement 

of the local-flow-inclination contours for both the predictions and the experimental 
data. 


CONCLaOING REMARKS 

In general it may be said that the two inviscid numerical methods investigated 
in this paper were capable of predicting forebody flows to levels consistent with the 
preliminary design of engine-inlet locations • Comparisons of contours of local 
angles of attack, local angles of sideslip, and local Mach numbers were compared with 
wind tunnel data. The supersonic code, STEIN, showed a marked sensitivity to any 
development of subsonic regions, especially at low supersonic Mach numbers or at 
angles of attack greater than 5®. The transonic code, WIBCO, was more robust and was 
able to calculate reasonable solutions through the Mach number range of 0.9 to 1.6 at 
angles of attack of 15® and over, depending on configurations. Both codes require a 
fairly elaborate geometric input, which may become time consuming. 

For aircraft inlet configurations more exotic than those considered, the two 
codes may not perform satisfactorily, particularly at high angles of attack or side- 
slip. Under such flight conditions, boundary layer build-up cannot be neglected even 
when the flow remains attached. This situation suggests that the next step in the 
analysis of such configurations should include viscous effects, including vortex 
shedding. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 2, 1984 
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( a ) Bow shock • 



(b) Canopy and wing shocks* 

Figure 2.- STEIN grid-shock development. 
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Figure 3.- STEIN cross-sectional grid pattern. 
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(a) a^. 


Prediction 

Experiment 



(b) 3^. 

Figure 7.- Contours predicted from STEIN 
code and Tailor-Mate experimental data 
for side-mounted configuration at 
= 2.5 and a = 0®. 
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(b) 3^. 

Figure 13.- Contours predicted from WIBCO 
code and Tailor-Mate experimental data 
for side-mounted configuration at 
































(b) 3^. 

Figure 16.- Contours predicted from WIBCO 
code and Tailor-Mate experimental data 
for side-mounted configuration at 























4 < 







Figure 18.- Contours predicted from WIBCO 
code and Tailor-Mate experimental data 
for side-mounted configuration at 
= 1 .6 and a = 

00 
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(b) 3^* 

Figure 19.- Contours predicted from WIBCO 
code and Tailor -Mate experimental data 
for side-mounted configuration at 
= 2.2 and a = 0®. 











(b) 3^. 

Figure 21 •- Contours predicted from STEIN 
code and Tailor-Mate experimental data 
for wing-shielded configuration at 
M_ = 2,5 and a = 0®* 







(c) 


Figure 21 






(b) 

Figure 22 •- Contours predicted from STEIN 
code and Tailor-Mate experimental data 
for wing-shielded configuration at 
M^ = 2.5 and a = 5°. 
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Figure 22 •- Concluded 
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